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ABSTRACT 

Substitution Delone set families are families of Delone sets X = {Xi, . . . , X^) which satisfy 
the inflation functional equation 

m 

Xi=\l {A{Xj) + Vij), l<i<m, 
i=i 

in which A is an expanding matrix, i.e. all of the eigenvalues of A fall outside the unit circle. 
Here the Pjj are finite sets of vectors in M.'^ and V denotes union that counts multiplicity. 

This paper characterizes families X = {Xi, Xn) that satisfy an inflation functional equa- 
tion, in which each Xi is a multiset (set with multiplicity) whose underlying set is discrete. 
It then studies the subclass of Delone set solutions, and gives necessary conditions on the 
coefficients of the inflation functional equation for such solutions X to exist. It relates Delone 
set solutions to a narrower subclass of solutions, called self-replicating multi-tiling sets, which 
arise as tiling sets for self-replicating multi-tilings. 
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1. Introduction 



Aperiodic and self-similar structures in M'^ have been extensively studied using tilings of 
as models. Among these are the classes of self-similar tilings and, more generally, self- 
affine tilings. Such tilings have been proposed as models for quasicrystalline structures (Q, 
Q, 0, [22 1, fS^j); they also arise in constructions of compactly supported wavelets and 



multiwavelets ( |l[,y], ||5[)- An alternate method of modelling quasicrystalline structures uses 
discrete sets, more specifically Delone sets (defined below), see |14], |15], |l^, [|^], [^], [28| 



|33]. These sets model the atomic positions occupied in the structure. In terms of tilings, such 
discrete sets can be viewed either as tiling sets, representing the translations used in forming 
tilings by translation of a finite number of different prototile types, or as control points marking 
in some way the location of each tile. 

Comparison of these two types of models, which appear rather different, motivates the 
question: Is there a notion of self-similarity appropriate to discrete sets and Delone sets? This 
paper develops such a notion, which is based on a system of functional equations dual to the 
functional equations associated to self-affine tilings and multi-tilings. 

We first recall the functional equation associated to the construction of finite sets of tiles 
{Ti, . . . ,T„} which tile R*^ with special kinds of self-affine tilings. The tiles are solutions to 
a finite system of set-valued functional equations which we call multi-tile equations, which 
encode a self-affine property. An infl,ation map (j) : M*^ —>■ R'^ is an expanding linear map 
(/)(x) = Ax in which A is an expanding n x n real matrix, i.e. all its eigenvalues |A| > 1. Let 
{Vij : 1 < i,j < n} he finite sets in R*^ called digit sets. 

Multi-Tile Functional Equation. The family of compact sets (Ti, T2, . . . , T„) satisfy the 
system of equations 

n 

A(r,) = |J(T,+p,,) (1.1) 



for 1 < i < n. 

This functional equation is set-valued, i.e. points are counted without multiplicity in the 
set union on the right side of (p. . 1|) . The subdivision matrix associated to ( |1.1[ ) is 

S=[|A,|]l<„-<n. (1.2) 

These functional equations have a nice solution theory when the substitution matrix satisfies 
the following extra condition. 

Definition 1.1. A nonnegative real matrix S is primitive if some power has strictly positive 
entries. 



In the one-dimensional case the subdivision matrix associated to ( |1.1| ) is always primitive. 
It is known that when S is primitive the functional equation (|1.1| ) has a unique nonempty 
solution 

T:= (ri,...,r„) 

in which all are compact sets (see |[l|, Theorem 2.3], [Q]); in this case all Tj are nonempty. 
In the imprimitive case there can be more than one nonempty solution T := (Ti, . . . ,r„) in 
which all Tj are compact sets (see examples in [^), and the theory becomes more complicated. 
There are however at most finitely many such solutions, and under general conditions there 
is a unique "maximal" compact solution, in which all Tj are nonempty, see Q. In this paper 
we primarily consider functional equations where the subdivison matrix S is primitive, though 
some results apply for general S. 

We are interested in the case where these sets Tj have positive Lebesgue measure and can 
be used in tiling M" by a self-affine tiling. The tiling sets for such tilings are special solutions 
to a second functional equation involving the same data, which is "adjoint" to the multi-tile 
functional equation. This functional equation counts multiplicities of sets, unlike (1.1), and we 
consider solutions to it that are multisets. 



Inflation Functional Equation. The multiset family X := (Xi, X2, . . . , X„) satisfies the 
system of equations 

n 

Xi=\l (A(X^) + Vij), 1 < i < n, (1.3) 
i=i 

where Dij are finite sets of vectors in M^. 

Here each Xi is a multiset (as defined in §2), and VJ=i denotes multiset union, as defined 
in §2. We can view the inflation functional equation as a fixed point equation = X, 

where ip{-) is an operator taking multisets to multisets, defined by ip{X) := {X^, ...,X^) with 

n 

X[:=\l{k{X,)+V,,). 

i=i 

The inflation functional equation can also be written as a system of equations for the multi- 
plicity functions: 

n n 

^xM) = E^A(X,)+A,(^) = E E ^x,i^-H^ - d)), for all x G M^. (1.4) 

i=i i=i deVij 
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A self-replicating multi-tiling consists of a pair (T, X) of solutions to the multi-tile functional 
equation and the inflation functional equation such that: 

(1) The solution T := (Ti, . . . ,r„) to the multi-tile functional equation has all sets Tj of 
positive Lebesgue measure. 

(2) The solution X = {Xi, . . . , Xn) to the inflation functional equation are sets (multi-sets 
with all multiplicities one) and Ur=i(^* + -^i) * tiling of M'^, using the T/s as prototiles. 

This notion extends the notion of self-replicating tilings studied in Kenyon , , which allow 
only one type of tile (n = 1). Later studies (|10|, |jll|) allowed n tile types but restricted the 



inflation matrix A to be a similarity. 

We call a family X a self -replicating multi-tiling family (for a fixed inflation functional 
equation) if it appears as part of a self-replicating multi-tiling {T,X). It follows that each 
multiset Xi in such a family is a uniformly discrete set. There are very strong restrictions 
on the data (A, Pjj) on the functional equations for a self-replicating multi-tiling to exist, for 
example the Perron eigenvalue condition given below. 

This paper studies solutions to the infiation functional equation that are discrete, including 
solutions that do not correspond to tilings. The inflation functional equation has properties 
that significantly differ from those of the multi-tile functional equation. For example, the multi- 
tile functional equation has a unique solution (in the case of primitive subdivision matrix) 
because the solution T is given by the unique fixed point of a contracting system of mappings. 
In contrast, the infiation functional equation involves an expanding system of mappings, and its 
solutions are not compact sets. It may have infinitely many different solutions, even infinitely 
many very nice solutions in some cases. Our replacement for the "contracting" condition, is to 
restrict attention to solutions X = (Xi, . . . ,X„) having the property that all Xj are discrete 
multisets. A set X in is discrete if each bounded set in M"^ contains finitely many elements 
of X. A multiset X is discrete if its underlying set X is discrete and each element in X has a 
finite multiplicity. A multiset family X is discrete if each multiset Xi in it is discrete. 

In §2 we give precise definitions and statements of the main results in the paper. Below we 
summarize the general contents of the other sections. 

In §3 we develop a structure theory for multiset solutions to the inflation functional equation 
that are discrete. These results are proved for general infiation functional equations, with no 
primitivity restriction. We show that discrete multiset solutions decompose uniquely into a 
finite number of irreducible discrete multisets, and that each irreducible discrete multiset is 
characterized by a finite set of points in it. However there are significant restrictions on the 
inflation functional equation data {l\,T>ij) necessarry for the existence of any discrete multiset 
solution, as indicated in later sections. 

In §4-§7 we study multiset solutions which correspond more closely to tilings. A multiset 
X is weakly uniformly discrete if there is a positive radius r and a finite constant m > 1 such 
that each ball of radius r contains at most m points of X, counting multiplicities; it is relatively 
dense if there is a radius R such that each ball of radius R contains at least one point oi X. A 
multiset X is a weak Delone set X if it is weakly uniformly discrete and relatively dense. 

A solution X = {Xi, Xn) to the inflation functional equation is a weak substitution 
Delone multiset family if each Xi is a weak Delone multiset; In studying solutions which 
are weak Delone multisets, we restrict attention to inflation functional equations that have a 
primitive subdivision matrix S. Then we make use of Perron- Probenius theory, which asserts 



3 



that a primitive nonnegative real matrix M has a positive real eigenvalue A(M) such that: 



(i) A(M) has multiplicity one. 

(ii) A(M) > |A'| for all eigenvalues A' of M with A^ A(M). 

(iii) M has both right and left eigenvectors for eigenvalue A(M) which have positive real entries. 

We call A(M) the Perron eigenvalue of M; it is equal to the spectral radius of M. 

For inflation functional equations whose substitution matrix is primitive, we show that 
those multiset equations having a weak Delone set solution must satisfy the following: 

Perron eigenvalue condition: The Perron eigenvalue A(S) of the subdivision matrix S 
satisfies 

A(S) = |det(A)|. (1.5) 

More generally, there are inequalities relating the Perron eigenvalue and properties of so- 
lutions of the multi-tile and inflation functional equations, as follows: 

(1) A necessary condition for the multi-tile functional equation with primitive subdivision 
matrix to have T = (Ti, ...,Tn) with some (and hence all) Tj of positive Lebesgue measure is 
that 

A(S) > |det(A)|. 

(2) A necessary condition for the inflation functional equation with primitive subdivision 
matrix to have a solution X = (Xi, Xn) with some (and hence all) Xi weakly uniformly 
discrete is that 

A(S) < |det(A)|. 

Inequality (1) is established by taking the Lebesgue measure on both sides of the multi- 
tile equation ( |l.lD , see Theorem 5^. In the special case when A is a similarity Mauldin and 



Williams pq] showed the stronger result that if A(S) < |det(A)| then the Hausdorff dimension 



of each Tj must be less than d. Inequality (2) is established in Theorem . 

In §5, under the assumption that the Perron eigenvalue condition holds, we give some 
necessary and sufficient conditions for the associated multi-tiling equation to have solutions 
Ti having positive Lebesgue measure. Then in §6 we use these to show that this condition 
is equivalent to existence of weak substitution Delone set solutions = X) is closely 

associated with self-replicating multi-tilings. In particular, the associated multi-tile functional 
equation necessarily has a compact solution T = (Ti, . . . , T„) with the Tj all having positive 
Lebesgue measure, see Theorem 2.4, and some iterate ip^{-) of the inflation functional equation 
then has a solution that is a self-replicating multi-tiling. In §7 we supplement this with a 
sufficient condition to have a self-replicating multi-tiling. 

In §8 we give examples showing the limits of our results. 

We distinguish multi-tilings, which are tilings of M"^ using several types of tiles, from multiple 
tilings, which are arrangements of tiles in M'^ such that almost all points in are covered 
exactly p times, for some positive integer p. We call p the thickness of the multiple tiling, and 
sometimes call such a tiling p-thick. A 1-thick tiling is just an ordinary tiling. The constructions 
of this paper can have associated to them tilings which are multiple in both senses, that is, 
p-thick multi-tilings for some p > 2. By definition self-replicating multi-tilings are 1-thick 
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tilings. We remark that the tiling set of a p-thick tiling cannot in general be partitioned to 
give a union of 1-thick tilings, even in dimension 1. 

To conclude this introduction, given any Delone set X one can associate a topological 
dynamical system ([[X]],R'^) with an R'^-action, in which [[X]] is the closure of the orbit of X 
under the M'^-translation action in an appropriate topology, see Solomyak |3^] . For substitution 
Delone set families these dynamical systems can be viewed as a generalization of substitution 
dynamical systems (Queffelec |3^), in that every primitive substitution dynamical system 
is topologically conjugate to a suitable substitution Delone set dynamical system. Under 
sufficiently strong extra hypotheses substitution Delone set dynamical systems are minimal 
and uniquely ergodic. We hope to return to this question elsewhere. 

Acknowledgments. We are indebted to the reviewers for very detailed comments improving 
the paper. 

Notation. The positive integer d refers to the dimension W^. We use -B_r(x) := {y G M"^ : 
||y — xjl < R} to denote the closed Euclidean ball of radius R around x. The positive integer 
n denotes the number of tile types or terms in the inflation functional equation. The positive 
integer p refers to the period of a cycle, and also sometimes to the the thickness of a multiple 
covering or tiling of W^. 

2. Definitions and Main Results 

We consider solutions to the inflation functional equation that are multisets. 

Definition 2.1. (i) A multiset X in is a set X_ together with a positive integer valued 
function, the multiplicity function mx '■ 2L ~^ ^>0) which assigns to each element x £ A 
a "multiplicity" mx (x) . Given a multiset A we shall use A to denote the underlying set 
of A, i.e. without counting the multiplicity, and use either mx(x) or m(A, x) to denote 
the multiplicity function. A multiset is called an ordinary set if every element in A has 
multiplicity 1. We extend the multiplicity function to all x € M'^ by setting mx(x) = 
if X ^ A. 

(ii) Given multisets Ai and A2 in we say Ai C A2 provided A^^ C A2 and mx^ (x) < 
mXalx) for all x G A2. In particular, A C A. 

It is convenient to extend the definition of multiplicity function mx(x) to all x £ M*^, by 
setting mx (x) = if x A . 

The term "multiset" is attributed to N. G. deBruijn (see Knuth [|^, p. 636]) and the 
concept has many uses. The combinatorial view of a multiset is as a collection of possibly 
repeated elements. For example, A = {0,0,1,3,4,4,4} represents a multiset in M in which 
is counted twice and 4 is counted 3 times. Thus A = {0, 1, 3, 4} and 

mx{0)=2, mx{l) = mx{3) = l, mx(4) = 3. 

One may alternatively regard a multiset A in M*^ as a nonnegative integer-valued pure 
discrete measure on W^. The results of this paper could be formulated in a measure-theoretic 
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framework, and the inflation functional equation then expresses an equality of measures. A 
reviewer observed that a measure-theoretic treatment might prove useful for further general- 
izations, but we do not attempt it here. 

Definition 2.2. For any multisets X and Y the multiset union XM Y is the multiset having 
the multiplicity function 

mxwY = mx + my, (2.1) 
and the multiset intersection X AY is the multiset having multiplicity function 

mxAY = min{mx,mY}. (2.2) 

For a multiset X and a set (or multiset) V the multiset sum X + V is 

X + V:= y (X + d). (2.3) 

Definition 2.3. A multiset family X = (ATi, . . . ,X„) is a finite vector of multisets Xi. Wc 
call a multiset family X = {Xi, . . . , X„) an n-multiset family. For n-multiset families X^^^ and 
we define 

x^^) V -^(2) = (xf ) V xf\. . . ,xW V 4^)). 

Definition 2.4. A multiset family X = {Xi, ...,X„) is discrete if for each i the multiset Xi is 
discrete, i.e. the underlying set 2Li is discrete and elements in Xi have finite multiplicity. 

In §3 we develop a structure theory for discrete multiset families, which decomposes them 
into irreducible families. Given an inflation functional equation, let be the inflation op- 
erator associated to it, which takes n-multiset families to n-multiset families, as defined in 
§3. 

Definition 2.5. (i) A multiset family X = {Xi, X„) satisfying an inflation functional equa- 
tion tp{X) = X is indecomposable if it cannot be partitioned as X = X^^^ W X^'^^ with each AfW 
nonempty such that ijj^X^^^) = X^^^ for each i = 1,2. Otherwise it is decomposable. 

(ii) A multiset family X = (Xi, X„) satisfying an inflation functional equation ip{X) = X 
is irreducible if there does not exist a nonempty X' = (X^, ...,X^) satisfying ip{X ) = X with 
X C X in the sense that X- C Xi for all i, with some X^ ^ Xi. Otherwise it is reducible. 

The concepts of indecomposable and irreducible multiset family (for a flxed inflation func- 
tional equation) are equivalent. If X is irreducible then it is necessarily indecomposable. 
Conversely, if X is reducible, with X' C X, and ip{X') = X' , then X" = X — X' (subtraction 
done on multiplicity functions) is a multiset family with X = X'y X", and linearity of the 
inflation functional equation yields ip{X") = X", so that X is decomposable. 

Theorem 2.1 (Decomposition Theorem) Let X be a multiset family X that is discrete 
and satisfies an inflation functional equation ip{X) = X. Then X uniquely partitions into a 
finite number of irreducible discrete multiset families that satisfy the same inflation functional 
equation. 
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An inflational functional equation il'i'^) = may not have a discrete solution, see Example 
8.4. In §3 we characterize indecomposable discrete multiset families X satisfying inflation 



functional equations, as follows. 

Theorem 2.2. Let X be a multiset family which satisfies an inflation functional equation 
= X and is discrete and indecomposable. Then X is irreducible and is generated by a 
finite "seed" S^^^ = 52, which consists of a periodic cycle Y = {(xj, ij) : 1 < j < p} 
in which xj G Xi. and there is some dj G 'Di-^-^^i. with Xj+i = Axj + dj, with (xp+i,ip+i) = 
(xi, ii). That is, 

X = hm 

The periodic cycle Y is the only periodic cycle in X and its elements each have multiplicity 
one. 



This result appears, in a more precise form, as Theorem 3^. For each p there are only 
finitely many periodic cycles Y and they can be effectively enumerated. However not all 
periodic cycles Y generate irreducible discrete multiset families. We show that there is an 
algorithmic procedure, which when given any such "seed" as input, has one of three outcomes: 

(1) If the generated multiset system X is discrete and irreducible, the procedure eventually 
halts and certifies this holds. 

(2) If the generated multiset system X is not irreducible, or if the limit does not exist, the 
procedure eventually halts and certifies this holds. 

(3) If the generated multiset system X is irreducible and not discrete, the procedure may 
not halt. 

We also prove in §3 a dichotomy concerning the multiplicities of elements appearing in 
the multisets Xi in an irreducible discrete multiset family X satisfying an inflation functional 
equation having a primitive subdivision matrix : Either they are all bounded by the period p 
of the generating cycle or else they all have unbounded multiplicities. (Theorem ^^.) We do 
not know of an effective computational procedure to tell in general which of these alternatives 
occurs. 

In §4-§7 we restrict to the case of inflation functional equations having primitive subdivision 
matrix. We study Delone set solutions and self-replicating multi-tilings. 

Definition 2.6. (i) A multiset X is weakly uniformly discrete if there is a radius r > and 
a finite m > 1 such that any open ball of radius r contains at most m points of X, counting 
points with multiplicity. 

(ii) A multiset X is uniformly discrete if there is a radius r > such that such that any 
open ball of radius r contains at most one point of X, counted with multiplicity. Such an X 
is necessarily an ordinary set. 

(iii) A multiset X is relatively dense if there is a finite radius R > such that each closed 
ball of radius R contains at least one point of X. 

Note that a finite union of weakly uniformly discrete sets is weakly uniformly discrete. 

Definition 2.7. (i) A multiset X is a weak Delone multiset if it is weakly uniformly discrete 
and relatively dense. 
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(ii) A multiset X is a Delone set if it is uniformly discrete and relatively dense. That is, it 
is an ordinary set which has both a finite nonzero packing radius and a finite covering radius 
by equal spheres. 

A finite union of Delone sets need not be a Delone set, while a finite union of weak Delone 
multisets is a weak Delone multiset. 

Definition 2.8. (i) A multiset family X = {Xi, . . . is a weak Delone multiset family if 
X satisfies the inflation functional equation (|l.3| ) in and each multiset Xi is a weak Delone 
multiset. 

(ii) It is a substitution Delone set family if it is a weak substitution Delone multiset family 
and each Xi is a Delone set. 

Weak substitution Delone multiset families can sometimes be associated to p-thick multiple 
tilings, for some p > 2. 

In §4 we prove the following result. 

Theorem 2.3 (Perron Eigenvalue Condition) If the inflation functional equation ipi^X) = 
X with primitive subdivision matrix S = [jl'tjl] has a solution X that is a weak Delone multiset 
family, then the Perron eigenvalue A(S) of S satisfies 

A(S) = |det(A)|. 

The narrowest class of solutions to the inflation functional equation are those having the 
following tiling property. 

Definition 2.9. A multiset family X = (Xi,...,X„) is called a self-replicating multi-tiling 
family for a given system (A, Pj^) in ( |I.3| ) with primitive subdivision matrix S if 

(i) is a substitution Delone set family for (A, Pjj). 

(ii) The associated multi-tile equation (0) has a unique solution T := (Ti,...,T„) with 
each Tj of positive Lebesgue measure. 

(iii) The sets {Tj + Xi : 1 < i < n and Xj G X^} tile W^. 

This definition requires that the multiset family X give a 1-thick tiling, rather than a 
p-thick multiple tiling for some p > 2, or a p-packing for some p > 1 that is not a tiling. 

In §5 we study the multi-tiling functional equation and show that a necessary condition for 
a solution with some Tj having positive Lebesgue measure is that A(S) > | det(A)|. Assuming 
that the Perron eigenvalue condition A(S) = | det(A)| holds, in Theorem |5.5| we characterize 
when solutions of positive Lebesgue measure exist. 

In §6 we prove the following result, which shows that substitution Delone set families are 
related to the existence of self-replicating multi-tilings. In this result {■) denotes the inflation 
operator composed with itself N times, whose associated inflation matrix is A''^. 
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Theorem 2.4. Let ip{X) = X he an inflation functional equation that has a primitive subdi- 
vision matrix. Then the following conditions are equivalent: 

(i) For some N > there exists a weak Delone multiset family X satsfying 'ip^{X) = X. 

(ii) For some N > there exists a self-replicating multi-tiling family X for the inflation 
functional equation if)^ {■). 

(iii) The subdivision matrix S satisfies the Perron eigenvalue condition A(S) = |det(A)| and 
the unique compact solution (Ti, . . . ,r„) of the associated multi-tile functional equation 
consists of sets Tj that have positive Lebesgue measure, I < i < n, and each Ti is the 
closure of its interior Ti = T° . 



These conditions imply that the Perron eigenvalue condition holds. If we instead assume 
that the Perron eigenvalue condition holds then condition (i) can be relaxed to assuming only 



the existence of a weakly uniformly discrete solution, as shown in Theorem 5.1. It is known 
there are many restrictions on the data {A,T>ij) in order for the conditions of Theorem 2A to 
hold. A number of different necessary (or sufficient) conditions on the data {A,T>ij) are given 



in Kenyon Lagarias and Wang H]- [||, Praggastis [||], H and Solomyak |§"|32 

In §7 we give a sufficient condition for a substitution Delone set family to be a self-replicating 
multi-tiling (Theorem |7.1| ). In §8 we give examples and counterexamples showing the limits of 
our results. 



3. Structure of Discrete Multiset Solutions 

We consider multiset families X = (Xi, . . . , Xn) satisfying the inflation functional equation 

n 

X,= \l iA{Xj) + Aj), l<i<n. (3.1) 

The inflation operator tp{-) associated to data {A,T>ij) is an operator which maps arbitrary 
n-multiset families X to n-multiset families ip{X) = X , as follows. Given X (not necessarily 
satisfying (|3.lD ), we define il^{X) = {X^, . . . , X„) by 



X'i := y {A{Xj) + Aj), 1 < i < n. (3.2) 



The inflation functional equation ( |3.lD then asserts that is a fixed point of V'(')) i-^- V'('^) = 
X. 

In this section we determine the structure of multiset families X that are discrete and satisfy 
ij^iX) = X, where we put no restriction on ip{-). The results allow imprimitive substitution 
matrices, unless otherwise stated. We also describe the structure of such sets that are also 
irreducible (Theorems 3.3 and 3.5). We show that every multiset family X that is discrete and 



satisfies ip{X) = X uniquely partitions into a finite number of irreducible sets (Theorem 2.1). 

We consider solutions to the inflation functional equation ( |3.l[) built up by an iterative 
process starting from a finite multiset which is a "seed." Let 

5W:=(5f\...,5f) (3.3) 
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be a system of finite multisets, and iteratively define the finite multiset family 5'^'^^^^ by 
^ i.e. 

n 

5f+i):=\/(A(5f) + P,,), l<i<n. 

We say that 5*^'^^ satisfies the inclusion property S^^^ C ^S*^"*^^ if 

g{0) ^ g{i) 1 < i < n , (3.4) 

in the sense of multisets. 

Lemma 3.1. // a multiset family S^^^ is finite and satisfies the inclusion condition S^^^ C 
?/;(5(o)) then 

^(fc) ^ ^(fc+i) ^ > Q _ ^3 5) 



The following limit set is well-defined 



Xi := lim Sf\ l<i<n , (3.6) 



fe^oo 



where X_^ is regarded as a countable point set. For each x G X ^ the multiplicity function has a 
limit, 

mxi(x) := lim m(5-''\x), (3.7) 

k~*oo 



where this multiplicity may take the value +oo. // all multiplicities in (3/j_) remain finite, then 



the multisets Xi are well-defined and the multiset family X = {Xi, . . . , X„) satisfies ^^{X) = X . 

Remark. If all the multiplicities remain finite, we write X := \\m.]^^^S^^\ The limit mul- 
tiset family X need not be discrete, and the underlying sets could even be dense. 

Proof. The inclusion (|3.5| ) follows by induction on k, since 

c yields 



The other statements follow easily from ( |3.5D . I 

To further analyze multiset families X that satisfy 'il^{X) = X, we put the structure of a 
colored directed graph on X. The vertices of this graph are the points in the disjoint union 
of the underlying sets Xj, with vertices being labelled (x, i), where x G 2Li ^ '^'^ and i is the 
"color". For each Xj € Xj and each dj G Dij, set 

X. := Axj + dj G Xi , 

and put a directed edge Xj — > xJ in the graph. We call x^ an offspring of Xj, and call xj a 
preimage of x.. We denote this (infinite) colored directed graph by G{X). We also assign to 
each vertex x G X_i & weight which is its multiplicity (x) . 

Lemma 3.2. Let X = (Xi,...,X„) be a multiset family which is discrete and satisfies the 
inflation functional equation ip{X) = X . 
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(i) There is a finite multiset family S^^^ Q X with the inclusion property, such that 

X = lim S^''^ . (3.8) 

fc— >oo 

(ii) The directed graph Q{X) contains a directed cycle. 

Proof, (i) We first show that there exists a radius R! > 0, depending only on such 
that every vertex in X can be reached by a directed path in the graph G{X) from some point 
xe X n Br,{0), where Br,{0) := {x € M'' : ||x|| < R'}. To show this, we set 

A(A-^) =niin{|A| : A is an eigenvalue of A (3.9) 

and A(A"^) < 1 since A is expanding. Fix a p > 1 such that A (A ^) < 1/p < 1. It is shown in 
Lind p3| that there exists a norm || • ||^ on with the property that 

P"^x||a < ^||x||a for all x € M'^ (3.10) 

The norm is defined by ||x||^ := X^fcli /c'^||^~^x||. 

We claim that there exists an i? > such that for any vertex x^ G with |[x^||^ > R, 
any preimage Xj G J£j of x^ must satisfy 



till A ^ ll^il 



To see this, let 



C = max{||dj||A : d^eVy, l<i,j<n} and := ^^C. (3.11) 



It follows from Xj = A (x^ — d^) that 

K-IIa<^(NIIa + c)<y^NIIa- (3-12) 

Thus if a vertex x has ||x||^ > R then every preimage of x in the graph G{X) is smaller in 
the norm |[ • ||^ by a multiplicative constant < 1. It follows that every vertex x in G{^) 
can be reached by a finite directed path in Q{X) starting from some vertex x' with ||x'j|^ < R. 
Since the set 

Br := {y gR^: ||y||A<i?}, 

is compact, it is contained in some ball of radius R' around 0, which gives the result. 

We now take the multiset family S^^^ = {S^\...,Sn^) such that Si consists of those 
elements x of Xi with ||x||^ < R, counting multiplicities, i.e. we take 

m (o)(x) := mxi(x) if x £ sf\ 

i 

The multiset family S^^^ is a finite family since X is discrete. The inclusion property 

5(0)C5W=V'(5W) (3.13) 
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holds because all preimages of vertices of Gi^^) in {x € M*^ : ||x||^ < R} lie in this set, since 
IIxjIIa <^{R + C) <R. Since C X and tp{X) = X we obtain 5^^) C X for all /t, viewed 
as multisets, so that 

X' = lim S'^''^ 

fc— >oo 

exists and is a multiset. Now S^^^ contains all points x of with the correct multiplicities for 
||x||^ < R. Now, by induction on /c > 0, one proves: 

For all X with ||x||^ < (^i±£) R, one has for 1 < i < n that the multiplicity 

m(5^''\x) = mx,(x). 

The base case A; = holds by construction, and the induction step follows because all 
preimages of a vertex in {||x|| < (i^)*^/?} lie in {||x|| < {■^-^)''~'^ R} . We conclude that 
X' = X, which gives (i). 

(ii) We choose R as in (i). Pick a point x in Af with 

xG^R:={yGM'^: [|yllA<i?}, 

which must exist by the argument in (i). All preimages of x necessarily lie in Bn. Since X 
is discrete, if we follow a chain of successive preimages of a vertex x in Q{X) we stay in the 
finite set X n Br, so some vertex must occur twice. The path from this vertex to itself forms 
a directed cycle in Q{X). I 

Recall that a multiset family X satisfying an inflation functional equation is indecomposable 
if it has no nontivial partition X = Xi V X2, with both Xi satisfying the same functional 
equation, and it is irreducible if there is no nonempty X' X satisfying the same functional 
equation. 

Theorem 3.3 (Irreducible Set Characterization) Suppose that the multiset family X is 
discrete and satisfies an inflation functional equation ip{X) = X . The following are equivalent: 

(i) X is indecomposable. 

(ii) X is irreducible. 

(iii) The graph Q{X) contains exactly one directed cycle (which may be a loop), and the 
elements of this cycle have multiplicity one. 

If condition (iii) holds, let Y = {xi,...,Xp} he the cycle in Q{X) with multiplicity one and 
define the multiset family 5^'^^ = {S^\ . . . , Sn^) by Si = {'x. ^ Y : x has color i}. Then 
C '0(5(°)) and 

X = lim S'^'^l (3.14) 

fc^oo 



Remark. In view of (iii) and (3.14) we call S^'^^ the generating cycle of the irreducible multiset 



system X, and we call p the period of this cycle. 
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Proof. We show (i) (ii) =^ (iii) =^ (i). Just before stating Theorem 2.1 we showed that 
imphcation (i) <^ (ii) holds for all multiset solutions to a fixed inflation functional equation. 

(ii) =^ (iii). Lemma |3.2| shows that X is generated by the points x in X that lie in the 
compact region ||x||^ < R, and that G{X) contains vertices forming a directed cycle Y inside 

this region. The vertices of Y define a finite multiset family 5^0) = (5(°\ . . . , Sn^), with sf'^ = 

n Y, and all elements of sf'^ have multiplicity one. It is clear that 5^°) C ^^^^ = ^(5(°)) 
since each element of S^^^ has a preimage in 5(0). Since each .S^^) C X the limit 

X' := lim 5(^) (3.15) 

k^oo 

exists and is a multiset family with X' X. We have X' = X , for if not, this contradicts the 
irreducibility of X. 

We next show that X has multiplicity one on the vertices of Y . Suppose not. The elements 
of 5(0) have multiplicity one, and at some later stage some vertex (x, i) of the cycle Y has 
multiplicity exceeding one in some S^'^\ Thus there exists a path to (x, i) of length k arising 
from some vertex (y, j) in 5(0), which does not arise purely from moving on the cycle Y . Since 
the vertex (y,j) can be reached from (x, i) by moving around the cycle Y, taking / steps, 
say, we obtain a directed path from (x, i) to itself of length k + / which does not stay in the 
cycle Y . Now we obtain two distinct directed paths of length (A; + l)p from (x, i) to itself, 
one by wrapping around the cycle Y k + I times, the other by repeating the path of length 
k + 1 p times. Since we can concatenate these two distinct paths in any order, it follows that 
the multiplicity of (x, i) in ^CC^+Opm) is 

at least 2™. This implies that the vertex (x, i) has 
unbounded multiplicity as m — > 00, so that lim^^oo "^^^^ does not exist, a contradiction. Thus 
X has multiplicity one on the cycle Y . 

Finally we show that X contains only one directed cycle. Suppose not, and that Q{X) 
contains a directed cycle Y' different from Y . The same argument as above says that X is 
generated by the cycle Y' and that all elements of Y' in X have multiplicity one. By exchanging 
Y and Y' if necessary we may suppose that Y has a vertex not contained in Y' . Now following 
a path from a vertex in Y\Y' to Y' we see that there exists an x' E Y' that has a preimage not 
in Y' . But x' also has a preimage in Y' because Y' is a cycle. This implies that the multiplicity 
of x' is at least 2, a contradiction. Thus Q{X) contains exactly one directed cycle. 

(iii) =^ (i). Suppose that G{X) contains exactly one directed cycle y, of multiplicity one. 
We argue by contradiction. Suppose that X were not indecomposable. Then we can write 
X = X'y X", where both X' and X" satisfy the inflation functional equation. Now Lemma |3.2| 
(ii) implies that both Q{X') and Q{X") contain a directed cycle. Therefore Q{X) either contains 
two directed cycles, or else contains a directed cycle of multiplicity at least two. This contradicts 
the hypothesis. I 

We can now prove the Decomposition Theorem. 



Proof of Decomposition Theorem 2A^. If X is not irreducible then it is decomposable 
and we have X = X'^^"' V Af(2) where each X'^^ satisfies X^^ = ip{X'^^). 
that X and <Y(') are generated by elements in the bounded region ||x||y 



Lemma 3^ showed 
< R. So the total 



multiplicity of elements of each ,%'(*) in the region ||x||^ < R is strictly smaller than that of 
elements of X in the region ||xj|^ < i?. If one of Af(*) is not irreducible then we can further 
partition it, and decrease the total multiplicity in ||x||^ < R again. Since X is discrete, this 



13 



process will end in finitely many steps, yielding 



N 

1=1 

To see that the partition is unique, to each element Xi of the partition is associated a unique 
directed cycle in of multiplicity one. But each directed cycle in Gi^V), counted with 

multiplicity one, must appear in some element of each partition. I 



Another consequence of Theorem |3.3| is a finiteness condition on discrete irreducible multiset 
familities. 



Lemma 3.4. Let '4>{-) be an inflation functional equation with data (A,D,j). Then for each 
p > I there are at most finitely many discrete multiset families (possibly none) X satisfying 
ip{X) = X that are irreducible and have a periodic cycle in Q{X) of minimal period p. 



Proof. Any periodic cycle {(xq, io); (xi, ii), . . . , (xp_i, ip_i)} constructed in Theorem 3^ has 
the form 

Xfc = Axfe_i + dfc, l<k <p 
with dfc € T>i^^i^_^ and (xp,ip) = (xo,io)- It follows that 

xo = A''xo + f]AP-Mfc. 

fc=i 

Thus 

xo = -(Af - A""'d/c^ , (3.16) 

where the matrix A^ — I is invertible since A is expanding. There are only a finite set of choices 
for the digits {di, . . . , dp}, so the number of choices for xq is finite. I 

We continue the study of discrete multiset families X with V'('^) = that are irreducible. 
We can give some further information on the multiplicities that occur in such X, under the 
extra hypothesis that the subdivision matrix is primitive. 

Theorem 3.5 (Multiplicity Dichotomy) Suppose that the multiset family X is discrete 
and satisfies the inflation functional equation "ip{X) = X with primitive subdivision matrix. If 
X = (Xi, X„) is indecomposable, then exactly one of the following cases holds: 

(i) Every element of each multiset Xi has multiplicity at most p, the period of the generating 
cycle in the directed graph Q{X). 

(ii) For each multiset Xi in X , the multiplicities mxi(x) ofxGXi are unbounded. 



Remarks. (1) For the general case of imprimitive S, one can prove a weaker multiplicity 
dichotomy; it applies only to those Xi such that i lies in the strongly connected component of 
the underlying graph of S that contains the generating cycle. 

(2) Both cases of the Multiplicity Dichotomy occur, as shown by examples. Example ^.5| 
gives an occurrence of case (i) in which all Xi have bounded multiplicity, with some multiplic- 
ities exceeding one. 
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Proof. The multiphcity of a vertex (x, i) in the graph Q{X) is m{x,i) := mXi(x). Fix a 
point (xo,io) in the generating cycle. We claim: For any point (x, i) ^ (xo,io) in the 
multiplicity m{x,i) is equal to the number of cycle-free directed paths from (xo,io) to (x, z). 
This is proved by induction on the length k of the longest cycle-free directed path from xq to 
X, which is bounded by (|3.12D . The base case /c = 1 is immediate. We have 

{x',i)eP(x,j) 

where P(x,i) is the set of preimages of (x, i) in Q{X), Assuming (x, i) / (xo,io) all preimages 
have shorter longest cycle-free path, and the induction step follows. 

Now assume that some vertex (xi,j) has multiplicity m{xi,j) = q > p- Without loss 
of generality we may assume that (xi,j) and (xo,io) have the same color j = i^, for if not 
by the primitivity hypothesis we can always find a descendant (x']^,io) of (xi,j), and clearly 
m(x']^, io) > m(xi, j)) > p. By the pigeonhole principle there exist two cycle-free directed paths 
from (xo,io) to (xi,io), whose lengths are L and L' respectively with L = L' (modp). Say 
V = L + sp for some s > 0. We can then create a second directed path of length V by first 
going around the cycle s times in the beginning and then following the path of length L. We 
show that these two different directed paths force the multiplicity of vertices to be unbounded. 

Since each directed edge can be labeled by elements in T>ij, 1 <i,j < m, we label the two 
directed paths by 

Pi = (di,...,dz.), P2 = (d;,...,d^). 

The fact that xq and xi have the same color i implies that di E Pji, d/, G Vn. and d'^ E I'j'j, 
d^ E Vik' for some j, A;, j', k' . Evaluating the two paths yields 

L L 

XI = A^xo + A^'^dj- = A^^o + A^'^^i- (^-l"^) 

So X]j=iA''^^dj = A-'^^dj. This means that following the two paths Vi and V2 from 

(xi,io) will lead to the same vertex (x2,io), which also has color zq- This process can be 
continued to obtain vertices (xfc,zo), k > 1, all of which have the same color iq. Note that 
there are at least 2^ distinct directed paths from (xo,io) to (xfc,zo) as we may concatenate Vi 
and V2 in any combination. These 2^^ directed paths remain distinct after removing the cycles 
in the initial segment. Hence m(xfc,zo) ^ 2^^. So the multiplicity of Xi^ is unbounded. Now 
the primitivity of the subdivision matrix implies that any vertex (x, io) has a descendant in 
each Xi. It follows that the multiplicity function is unbounded on every Xi. I 

As we have shown, all irreducible discrete multiset families X are generated iteratively 
from a cycle as "seed" . However, the multiset family generated from a given cycle need not be 
discrete. The following lemma gives a criterion for discreteness. 

Lemma 3.6. Let be an inflation functional equation with data (A,Djj), and let S^^^ be a 
finite multiset family with the inclusion property S^^^ C S^^^ = il^{S^^^). Set R = R{A,'D) equal 
to the constant in l \3.1]\ ). If for some k > 1, 

S^''^ n {x E : ||x||a <R} = 5^'^+^) n {x E M'^ : ||x||a < R} (3.18) 

counting multiplicities, then the limit 

X := lim S^^^ 

k—*oo 
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is a discrete multiset family that satisfies = ■ Conversely, if X is discrete then (\3.l3i ) 

holds for all sufficiently large k. 



Proof. The iterative scheme ^S'-'^^"'^-* = tlj{S^^^) is geometricahy expanding outside the set 
{x G : ||x||a < R}, i.e. gives 

Ax + d > o x - C > X . 

p+1 

The condition ( |3.18| ) says that S^''^ stabihzes inside {x E R'^ : ||x||^ < R}. By induction on 
j >0 one obtains that S^'^'^^^ stabihzes inside the domain j^x € M"^ : ||x||^ < ^ '^^^^ ^ R^. The 
hmit family X is discrete inside each of these domains, hence is discrete. I 



Lemma S.C leads to a recursively enumerable procedure which recognizes all irreducible 
discrete sets X. More precisely, consider the class of inflation functional equations ^{■) whose 
data (A, Pjj) is drawn from a computable subfield K of C, which we take to be one in which 
addition and multiplication are computable, and one can effectively test equality or inequality 
of field elements. Then all finite cycles Y for such ip{-) will have elements in K. One can 
give a procedure which tests whether a given seed 5(0) generates an irreducible discrete X, 



using the criterion of Lemma 3.6. It will eventually halt for all such sets, and certify that they 
are discrete. (One does not need to compute i? or || • ||^ exactly to guarantee halting.) This 
procedure also detects all non- irreducible sets, by either finding a second cycle or a cycle with 
multiplicity. However this procedure is not recursive, for given an input cycle Y that generates 
an irreducible X which is not discrete, the procedure can run forever without halting. 

The following lemma gives a simple case when all irreducible solutions to the inflation 
functional equation are discrete (if any exist). 

Lemma 3.7. Suppose that the inflation functional equation tjj[-) has data (A, T>ij), in which A 
is an expanding integer matrix and all vectors in T>ij have rational entries. Then all irreducible 
multiset families generated by a periodic cycle are discrete. 



Proof. All cycles Y of period p are generated by solutions of the form ( 3.16| ), and the 



hypothesis guarantee that such a solution x is a rational vector. Starting from (x, i) one 
generates the orbit Y and a seed S^^^ which consists of rational vectors. All elements of S^^^ 
are rational vectors with denominators dividing dD, where d is the denominator of x and D 
is the greatest common denominator of all the rationals in Vij. Thus there are only a finite 
number of possible choices for rational vectors in the ball of radius R = R{A,'D) in ( 3.11| ). 



Thus the criterion of Lemma 3^ applies to conclude that if X = lim^^oo 5'-'^) exists, then it is 
necessarily discrete. I 



Lemma p.7| does not assert existence of irreducible multiset families, only that they are dis- 



crete if they do exist. Example 8.4 in §8 shows there are cases satisfying the hypotheses where 
no discrete solutions exist. At the other extreme, there are inflation functional equations of this 
type having infinitely many different irreducible discrete multiset families, see Example ^.3| . 

We conclude this section by raising two related (unsolved) computational problems: 

Computational Problem (i). Given a seed Sq which generates an irreducible discrete multiset 
family X satisfying a given inflation functional equation ip{X) = X and an index i, determine 
whether the multiplicities of points in Xi are bounded or not. 
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Computational Problem (ii). Given the same data as above, determine the maximum multi- 
phcity of points in Xi. 



4. Perron Eigenvalue Condition 

In this section we consider inflation functional equations having a primitive subdivision matrix 
S so that the Perron eigenvalue A(S) is defined. Our main object in this section is to show that 
a necessary condition for an inflation functional equation to have a solution that is a Delone 
set family is that the Perron eigenvalue S satisfies the Perron eigenvalue condition 

A(S) = |det(A)|, (4.1) 



as given in Theorem 2^. To achieve this we prove several preliminary results. Some of them 
apply more generally to multiset families X which are weakly uniformly discrete in the following 
sense. 

Definition 4.1. A multiset family X is weakly uniformly discrete if the union \/^^iXi is a 
weakly uniformly discrete multiset. 

A finite union of uniformly discrete multisets is weakly uniformly discrete if and only if 
each multiset itself is weakly uniformly discrete. So a multiset family X is weakly uniformly 
discrete if and only if each Xi in it is weakly uniformly discrete. 

We first consider an irreducible multiset family X satisfying X = ^^{X). Let S'^^^ be the 
generating cycle for X and set 5(^+1) = ^1^{S^^'^). Write S^^^ = 5^ \ . . . , 5^''^) and let 

(fc) Ik) 

' denote the total multiplicity of the multiset , which is given by 



Lemma 4.1. Let X be an irreducible discrete multiset family satisfying ^{X) = X for the 
data (A, Aj), with subdivision matrix S. Let H('=) = [i^f \ . . . , H^^Y ■ Then H('=) = S''H(o). 

Proof. We prove that h('=+^) = SH^'^). Observe that 

n 

Hence 



i=i 



yielding hC^+i) = SH'^'^l Therefore H^'^') = S^ll^°\ 
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Theorem 4.2. Let X he an irreducible discrete multiset family satisfying ip{X) = X for the 
data {A,'Dij), with primitive subdivision matrix S = Let -Bi(O) be the unit ball and set 



xGX,nA''(Bi(o)) 

Then there exist constants Ci,C2 > such that, for each k>l, 

CiA(S)'= < Mf ) < C2A(S)^ 



(4.2) 



(4.3) 



Proof. We first establish the lower bound. Observe that for each element x G VILi ^i'^'^ there 
exist digits di, . . . ,dfc in the collection of digit sets {Vij} such that x = o • • • o rdj.(xo), 
where xq is in the generating cycle and Td(x) := Ax + d. So 

k k 

X = A'^xo + ^ A^-M, = A^ (xo + ^(A-i)^+i-M,) . 
i=i i=i 

Since A~^ is contracting, there exists a constant C depending only on {Vij}, A and S^^^ such 
that 



xo + 5;(A-i)^+i-M, 



< C. 



Fix an integer fci > so that Bc{0) C A''i(5i(0)). Then x G A''+'^i(Si(0)). It follows that 
for any x G X ^ n k^^'^^ {Bi{Q)) we always have 

mx,(x) = m^{fc)(x). 

Therefore Mf''^^'^^ > H^^^ for all A; > 1. By Lemma 44 and the primitivity of S we have 
fC') > C[X{S)'' for some constant C[. The lower bound M^^') > CiA(S)'' follows by taking 
Ci = C[X{S)-^K 

To establish the upper bound we claim that there exists a k2 such that for each k > 1, each 
i and x G X ^ n A^(i?i(0)) and any vertex x G X ^ with x G A'^(i?i(0)) we must have 



mx,(x) = m^(fc+fc2)(x)- 



(4.4) 



The claim gives the upper bound, for (4^) implies that M-'^'' < Hf^^'^'^' holds for each k >1 
and each i. However Lemma 4.1 and the primitivity of S together show that there exists a 
constant Cn such that 



Ak+ki) 



Combining these inequalities gives the upper bound < C2A(S)^ on taking C2 = C^A(S)'=2. 



It remains to prove the claim (4^). We note that any x G 2Li ^ A^^Bi) must be in sf ^ for 



some / as a result of Lemma 3.1. Suppose that / > k. Then there exist some digits di, d2, 



,d, 



in Pjj's such that 



x = A^xo + ^A^-Mj, 
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and 

k 

xo = A-'x-^A-Mfc+i_,. (4.5) 
i=i 

Observe that xq is bounded since / > k, A^'^x E Bi and A is expanding. Hence there exists a 
constant i?o independent of k and / such that ||xo|| < Rq- Let 

n 

k2 = ^\x,nBR,{0)\. 

i=l 

We now argue by contradiction, and assume that (^^) is false. Then there exist some 
k > 1, i and x € X j n A^{Bi) such that mx,(x) > m (fc+fe2)(x). This means there exists a 

cycle- free directed path in G{^) of length / > k + k2 from an element xq in the generating cycle 
S^^^ to X, say x = tj, o • • • o ° ''"di(xo)- Let Xj = rj^. (xj_i), 1 < j < Z, be the vertices of 
this directed path. It follows that ||xj|| < Rq for < j < I — k. But there are only k2 vertices 
in i?/{(,(0) and I — k > k2, hence there exist two identical vertices among {xj : 1 < j < I — k}. 
This contradicts the assumption that the directed path is cycle- free, proving ( |4.4| ). I 



Theorem 4.3. Let X be a weakly uniformly discrete n-multiset family satisfying ip{X) = X 
for the data {A,Vij), with primitive subdivision matrix S.Then 

A(S) < |det(A)|. 



Proof. We first assume that X is irreducible. The weak uniform discreteness of X implies 
that M^''^ < CVol(A^5i(0)) = C|det(A)|'= for some positive constant C. From Theorem U 
it immediately follows that |det(A)| > A(S). 

If X is reducible, then X = VjLi '^^-''^ ■ The argument above now applies to X^^^ to yield 
A(S) < |det(A)|. I 



Proof of Theorem |2.3| . Since X is a weak substitution Delone multiset family, it is weakly 
uniformly discrete. By Theorem |4.3| this immediately gives A(S) < | det(A)|. 

To prove the other direction A(S) > | det(A)|, let X = V^i X^^\ Let mJ^^ be as in Q, 
but defined for X^^\ Then the relative denseness of each Xi yields 



N 



Y^^i j ^ C'Vol(A^'Bi(0)) = C7'|det(A)|^' 
i=i 

for some positive constant C. Hence max,- Mj;j' > ^C'\ det(A)|'''. Taking the A;-th roots and 
letting A; — > cxo now yields A(S) > | det(A)|, using Theorem |4.2|. I 



5. Multi-Tiling Functional Equation 

Our object in the next section is to relate the existence of weak substitution Delone multiset 
families to self-replicating tilings. In this section we obtain preliminary information, concerning 
solutions of the multi-tiling functional equation. Theorem f>.2\ gives a necessary condition for 
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solutions of positive Lebesgue measure and Theorem |5.5| gives a set of equivalent conditions for 
positive Lebsegue measure that apply when the Perron eigenvalue condition A(S) = |det(A)| 
holds. 

We begin with a basic existence result on solutions to the multi-tiling equation. 
Proposition 5.1. The multi-tile functional equation 

n 

A{Ti) = (j {Tj + Vji), 1 < i < n, (5.1) 

with primitive subdivision matrix S has a unique nonempty solution (Ti, . . . , T„) in which each 
Ti is compact. In this solution all Ti are nonempty, and 

oo 

Ti = { A-Hu,_, I d,,,,_, G (io, Ji, j2, . . .) G {1,2, . . . ,iVf ^ jo = i}- (5.2) 

k=l 

Proof. Flaherty and Wang Proposition 2.3] prove under the hypothesis 

n 

(*) U ■^■J* nonempty for 1 < i < 



that the multi-tile equation ( |5.1| ) has a unique nonempty solution (Ti, . . . ,T„) in which all Tj 
are compact sets, and that in this solution all Tj are nonempty. The primitivity assumption 
on S implies that hypothesis (*) holds, hence this result applies. I 

We remark that ( ^.2[ ) has a graph-theoretic interpretation. Form a directed graph Qi^) 
whose vertices correspond to sets and with a directed edge from T, to Tj labelled d if d G Vji. 
Then (^]^) is equivalent to saying that for each point x G Tj there exists an infinite directed 
path (di, d2, ds, . . .) in the graph Q{iIj) where an edge connects with di G Vji for some j such 
that 

oo 

x = ^A-'=dfc (5.3) 

k=l 

and vice versa. 

We next give a simple necessary condition for the sets Tj to have positive Lebesgue measure. 
Theorem 5.2. Suppose that the multi-tile functional equation 

n 

k{Ti) = y {Tj + Vji), l<i<n, (5.4) 

with primitive subdivision matrix S = has a nonempty solution (Ti, ...,Tn) consisting of 

compact sets, in which at least one Ti has positive Lebesgue measure. Then 

A(S) > |det(A)|. (5.5) 
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Proof. Taking the Lebesgue measure of both sides of the multi-tile equation ( |5.4| ) gives, for 
1 < i < n, 

n 

\detiAMT,)<Y,\^jMT,), 
i=i 

where {Vjil denotes the cardinality of Vji (counting multiplicity). By hypothesis v := [//(Ti), . . . , //(T„)] 
is a nonnegative row vector, not identially zero, and the equation above gives 

vS > |det(A)|v. 

Write vS = |det(A)|v + wi, where wi is nonnegative. Repeated multiplication by S and 
back-substitution yields vS'^ = | det(A)|^'v + w^,, where w/; is nonnegative. We conclude that 
the spectral radius p{S) satisfies 

p{S) > lim llvS^II^/'^ > |det(A)j. 

fe— +00 

Thus A(S) > |det(A)|. I 

In the remainder of this section we develop criteria for positive Lebesgue measures of the 
sets Ti that apply when the Perron eigenvalue condition A(S) = | det(A)| holds. 

We define the digit multisets 

n 

^^T^■■= V (^nr.-.+^^Jr.-^,m-.+■■■+^"'-'^.^), (5.6) 

jl,--;jm-l = l 

in which the sum is interpreted as counting multiplicities. It is easy to check that iterating 
(P) yields 

A-(r,) = U(T,+P-), i = l,...,n. (5.7) 
In (|5.7| ) we do not count multiplicity, so it suffices to use 'DJl instead of 2?^^. 

Definition 5.1. (i) A family of discrete multisets {£a '■ a G /} in M'^ is equi-uniformly discrete 
if there exists an eq > such that each Sq, is uniformly discrete and any two distinct elements 
in £a are separated by at least distance Eq. In particular, each is an ordinary set. 

(ii) A family of discrete multisets {£a ■ a € /} is called weakly equi-uniformly discrete if 
there exists a fixed M > such that for each a & I and each ball B of radius 1 in the 
number of elements of £a in B (counting multiplicity) is bounded by M. 



The following theorem is an extension of a theorem in Sirvent and Wang |34]. Its hypotheses 
can only be satisfied when the Perron eigenvalue condition holds, because the weak equi- 
uniformly discrete hypothesis can only hold when A(S) < | det(A)| 

Theorem 5.3. Suppose that the compact sets (Ti, . . . ,Tn) satisfy the multi-tile functional 
equation 

n 

A{Ti)= \J{Tj+Vji), i = l,...,n, 
i=i 

with primitive subdivision matrix S, and that at least one Ti has positive Lebesgue measure. If 
the collection of multisets {'^1 '■ 1 < i < ?ti > 1} are weakly equi-uniformly discrete, then 
each Ti has nonempty interior T?, and is the closure of its interior, Ti = T° . 
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The primitivity of the subdivision matrix impUes that if a single Tj has positive Lebesgue 
measure, then they all do. The proof of Theorem is based on the following covering lemma. 



Lemma 5.4. Suppose that (Ti, ...,Tm) are compact sets satisfying an inflation functional equa- 
tion with primitive subdvision matrix, and that at least one Ti has positive Lebesgue measure. 
Suppose also that the collection of multisets {^ji : 1 < i,j < n,m > 1} is weakly equi- 
uniformly discrete. Then given any sequence of positive numbers {6m} with limm^oo = 0, 
there exist positive constants Rq and Kq such that the following holds: For each m > 1 there 
exist subsets . . . ofM!^ contained in the ball BRg{0), each of cardinality bounded by Kq, 
such that := Uj=i(^j + ^J^) ^^-^ ihe property that 

fi{Bi{0) n Qm) > (1 - 5''+^6m)KBim- (5.8) 



Proof. Without loss of generality we assume Ti has positive Lebesgue measure, so Ti has a 
Lebesgue point x*, i.e. there is a sequence rm such that 

/i(S,„(x*) n Ti) > (1 - 5m)KBrA^*))- 

It follows that 

//f A'(B,„(x*) n Ti)) > (1 - <5^)/i(A'(B,„(x*))), for all / > 0. (5.9) 



We first show that for sufficiently large I = l{m), there exists a unit ball i?i(y) C A' (i?,.^ (x* ) ) 
with 

jBiiy) n A'(ri)) > (1 - 5''+^6m)KBim. (5.10) 



Indeed, since A is expanding, A\Br^{x*)) is an ellipsoid Oi^m whose shortest axis goes to 
infinity as / goes to infinity. Let O'l ^ be the homothetically shrunk ellipsoid with shortest axis 
decreased in length by 2, so that all points in it are at distance at least 1 from the boundary 
of Oi^rn- By a standard covering lemma (cf. Stein |^, p. 9]) applied to O^'^ there is a set 
{Si(y')} of disjoint unit balls with centers in O'l ^ that cover volume at least b~'^^{0[ ^). Also 
5~'^fi{Oi ,^) > 5~'^~^ fi{Oi^m) once the shortest axis is of length at least 2{d+ 1). All these balls 
lie inside Oi^m- By ( ^^ ) at most 5m/i(A'(i?r^(x*))) of the volume of A\Br^{x*)) is uncovered 
by A'(i3r„(x*) n Ti), so at least one of the disjoint balls {Bi{y')} must satisfy ( ^.10 ). 

By (^) we can rewrite the inequality ( 5.10| ) as 

n 

(i?i(y) n +2^i))) > (1 - 5'^'dmMB^m, 



whence 



i?i(o)n(U(r,+p^i-y)) I >(i-5'^+iUMi?i(o)). 



This shows that if we choose 

£-]":= {d - y I d € with (Tj + d - y) n 5i(0) / 0} 
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then (5.8) holds. Since all Tj are compact, all £j^ lie inside the ball BJi^^{0) with Rq := 
1 + maxj diam (Tj). The cardinality of all fj" are upper bounded by a constant Kq because the 
set of all is weakly equi-uniformly discrete. I 



Proof of Theorem 5.3. We apply Lemma 5.4 and choose a subsequence rrik — > oo so that 



|£'mfc| converges (as a multiset) for all j, and we denote the limit multisets by £j^. This can 
always be done because {f™} are uniformly bounded and have uniformly bounded cardinality. 
Clearly £j° has cardinality at most Kq. So 

n n 

^(i?i(o)n(U(T,- + £■-))) > iiminf^(si(o)n(|J(T, + £:7'=) 

j=i i=i 



> liminf(l-5'^+^5™J//(Bi(0)) 

fc^oo 



Since each Tj + is a closed set, we must have 

n 

B,{0)n[\J{Tj+£°^))=B,{0)- 
i=i 

This means that at least one of the Tj 's must have nonempty interior Tj . But if so then the 
primitivity of the subdivision matrix implies that all Tj have nonempty interior. Let Tj = Tj. 
Then (T(,...,T^) must also satisfy the same multi-tile equation (|5.1|). By the uniqueness 
assertion in Proposition we have Tj = Tj for 1 < j < n. I 



We obtain the following characterization of the positive measure property /u(T) > 0, which 
applies when the Perron eigenvalue condition holds. 

Theorem 5.5. Assume that the family of compact sets (Ti, . . . , T„) satisfies a multi-tile func- 
tional equation with the data (A, Pjj) with primitive subdivision matrix S satisfying the Perron 
eigenvalue condition A(S) = | det(A)|. Then the following conditions are equivalent: 



(i 



1111 
(iv 

(v 
(vi 

(vii 



For m > 1 all multisets T>Jl for 1 < i,j < n are ordinary sets, and the family {"DJ^ '■ 1 < 
i, j 1^ n, m > 1} is equi-uniformly discrete. 

For some fixed 1 < k,l < n and m > 1 the multisets are ordinary sets and the family 
{P^ : m > 1} is equi-uniformly discrete. 

The family of multisets : 1 < i, j < n, m > 1} is weakly equi-uniformly discrete. 

For some fixed 1 < k,l < n the family of multisets {P^^ : m > 1} is weakly equi-uniformly 
discrete. 

One set Tj has fJ,{Tj) > 0. 

Every set Tj has Tj = T° , hence all fi{Tj) > 0. 

Each Tj = T° , for 1 < j < n, and /i(5Tfe) = 0. 
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Proof. The implications (i) =^=- (ii) =^ (iv) and (i) =^ (iii) =^ (iv) are obvious. 

(iv) =^ (v). First it is easy to check that = [l^'jyl] where {V^Jl denotes the cardinahty 
(counting multiphcity) of P™. The primitivity hypothesis imphes that there exists a cq > 

such that, for all m > 1, \V^\ > coA(S)'" for ah 1 <i,j <n. Now for all i let 1]^°^ = Bi{0) 
and 

n 

i;M = UA-1(7;('"-^)+P,,), m>0. 

Then T^"^^ — > Tj in the Hausdorff metric (see [^). We prove that /a(T^) > for all 1 < / < n. 
To see this, we note that 

n 

A™(i;('")) = (J (Tj°) + p^^). 

So 

n 

A™ (t/™)) = U (Tf ) + P^"?) D rf ) + P^. (5.11) 

Since {T^^l} are weakly equi- uniformly discrete, there exists a constant ci > such that 
l^fcil ^ for all m. Now since are weakly equi-uniformly discrete there exists an 

M > such that each unit ball contains at most M elements of P^. Therefore each point in 
7^°^ can be covered by no more than M copies of T^^^ +d, d G P^. Therefore by (|5Tll ) 



|det(A)r/.(T/"^)) > ^\VTMTf'^) > ^5coCiA(S)-, 
where (5 is the volume of the ball -Bi(O). Using the hypothesis A(S) = | det(A)| yields 

M(r/'"^)>^<5coCi>0. 



It follows that, for all / 



n{Ti) > liminf^(r/"'^) > 0. 



(v) =^ (i). First we note that ^{Tj) > for all j as a result of the primitivity of S. Let 
e = [/i(Ti), /i(T2), . . . , /i(T„)]. Taking the Lebesgue measure of both sides of the iterated multi- 
tile equation (|5.7D yields 



n n 



A-(S)MT.)<^M7;)l2^l<2^Mr,)|P,"?|, l<i<n. (5.12) 

In other words, A™(S)e < eS™. But S™ is a primitive nonnegative matrix with Perron- 
Frobenius eigenvalue A'"(S). So ( 5.12| ) can hold only when e is a left Perron- Frobenius eigen- 



vector of S and all inequalities in ( |5.12D are equalities. This immediately yields \'DJl\ = 
hence each VJl is an ordinary set. The equalities in ( 5.12| ) also imply that the unions Tj 



are all measure-wise disjoint. So all P^, and hence all P^, are equi-uniformly discrete for 
some e > 0. 



(v) (vi). This follows from Theorem 5.3, since (v) also implies (i) 
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(vi) =^ (vii). We only need to prove that fj,{dTj) = for all j. Let v = [fi{dTi), . . . ,iJ,{dTn)]. 
We have for all i 



n n 



-jlJ 



Similar to (|5.12| ), taking the Lebesgue measure yields A'"(S)v < vS™. Again, this can occur 
only when v = or v is a Perron- Frobenius left eigenvector of S™". Assume that v 7^ then 
ah i^idTj) > 0, and 

n 

/.(9(A-(r,))) = l^T^l^'idTj). (5.13) 
i=i 

But Tj has nonempty interior, so for sufficiently large m > the inflated set A™'(Tj) will contain 
a sufficiently large ball in its interior. Since A'"(Tj) is the union of Tj +^™, 1 < i < there 
must be some k and d E P^^j such that T^ + d is completely contained in the interior of A™(Tj). 
Hence 

n 

5(A-(r,))c|J(5T,+p-)\(r, + d). 

So ( ^.13 ) is impossible, a contradiction. 

(vii) =^ (vi). This is obvious. I 



6. Weakly Uniformly Discrete Multiset Families and Tilings 

In this section we suppose that the Perron eigenvalue condition A(S) = |det(A)| holds. Using 
the results of §5 we relate the existence of positive Lebesgue measure tiles Tj for the associated 
multi-tiling equation to weakly uniformly discrete multiset solutions and self-replicating multi- 
tiling solutions to some iterate V'^(') of the inflation functional equation. Then we deduce 



Theorem 2.4. 



Theorem 6.1. Let = X be an inflation functional equation that has a primitive subdi- 

vision matrix S that satisfies the Perron eigenvalue condition 

A(S) = |det(A)|. 

Then the following conditions are equivalent: 

(i) For some N > there exists a weakly uniformly discrete multiset family X such that 
il>^{X) = X. 

(ii) For some N > there exists a uniformly discrete multiset family X such that ip'^ {X) = 
X. 

(iii) For some N > there exists a self-replicating multi-tiling family X with ip'^ (X) = X. 

(iv) The unique compact solution (Ti, . . . , T^) of the associated multi-tile functional equation 
consists of sets Ti that have positive Lebesgue measure, 1 < i < n. 
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Proof. Iterating N times the inflation functional equation X = ^^{X) on multisets gives a 
new inflation functional equation X = ip^ {?(!), which corresponds to 

n 

X,= \/{A''{X,)+Vf^), l<i<n, (6.1) 
i=i 

where the sums are interpreted as multiset sums. 

We show that (i) =^ (iv) =^ (iii) =^ (ii) =^ (i). To begin, the implications (iii) =^ (ii) =^ (i) 
are obvious. 

(i) =^ (iv). By assumption there exists an M > such that any unit ball in contains at 
most M elements (counting multiplicity) of each Xi. Replacing by mN in ( |6.lD it follows 
that any unit ball contains at most M elements of each T>^^ . Observe that by ( p.6p we have 

n 

1=1 

Therefore any unit ball contains at most M elements of each T^^j'^ if it contains at most M 
elements of each Dfj. This immediately yields the weakly equi- uniform discreteness of all the 
sets {V'^j}. So by Theorem ^ 

(iv) =^ (iii). Since T° $ it follows from the observation (|5.3| ) that there exists an infinite 
directed path (di,d2,d3, . . .) in the graph G{X) with di € Vji for some j such that 

oo 

xo = J^A-^dfc GTf. 

k=l 

Since all Dij are bounded, there exists an A^' > such that for all infinite directed paths 
(d'^, dg, dg, . . .) with d'j = dj for j < N' we also have 

oo 
k=l 

The primitivity of the subdivision matrix S now implies that we can find an infinite directed 
paths (d^,d2,d3, . . .) which has d* = dj for j < N' , and which is periodic for some period 
iV > TV' in the sense that d^_^^ = for ah k. Let xj^ = EfcLiA^'^d^,. Then x^ G Tf and 
A^xS = xS + d for d = J2fjo^ AM^_^.. Observe that d G P^. So we have 

-xSeA^(-xS) + 2?fi. (6.2) 

Now consider the inflation functional equation X = xp^{X). Set X^^^ = {xf'\ . . . ,Xn ^), 
with 

xf) = {-xS}, = 0, 4°) = 0. 

Define = V^(;f(™-i)) = V™^('^(°)). Set 

n 

X« = V(A"^(^f ) = A^({xS}) 
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hence we have x[^^ C x[^^ by ( |6.2D , and obviously we have = ^ C xf ^ for i > 2. So the 
inclusion property ;f(°) C V^(<^) = -^^^^ holds. It follows that C X^^^ C ^"(2) c .... 
But notice that 



We conclude that each X^"^^ is an ordinary set and is eo-uniformly discrete. Let Xi = 
U^=o andX = {Xi,..., Xn). Then X = and X is eo-uniformly discrete. 

It remains to show that A" is a Delone family and to establish the tiling property of X. 
Observe that is in the interior of Ti — Xg . Now 



r^(ri - X*) = U (T, + - A"^x*) = y (r, + x^""^] 

j=i i=i 



Note that the union IJj=i(^j + ^ji^) is measure-wise disjoint as proved in the proof of (v) 
=^ (i) in Theorem |5.5|. Taking the limit as m ^ oo will keep the union measure-wise disjoint. 



Hence Uj'=i(^j + Xj) is a tiling of R . So all Xj must be relatively dense as a result of the 
primitivity of S. This completes the proof. I 



We conclude this section by deducing Theorem 2.4 as a consequence of Theorem 6.1 



Proof of Theorem |2.4| . (iii) =^ (ii). This follows from Theorem |6.1| (iv) =^ (iii). 

(ii) =^ (i). A self-replicating multi-tiling family is a weak substitution Delone multiset family. 

(i)=^ (iii). By Theorem p.3| the inflation functional equation satisfies the Perron eigenvalue 
condition. A weak substitution Delone multiset family is a weakly uniformly discrete multiset 
family, and satisfies the Perron eigenvalue condition by Theorem p.3| . Thus the conditions of 
Theorem |6.1| (i) hold, and the result then follows from Theorem |6.1| (iv). I 



7. Self- Replicating Multi-Tiling Families 



In this section we study self-replicating multi-tiling families as a subclass of weak substitution 
Delone multiset families. 



Theorem 7.1. Let X he an irreducible weak Delone multiset family satisfying the inflation 
functional equation V'('^) = X for the data (A, T>ij), where the subdivision matrix S is primitive. 
Suppose that the fundamental cycle of X has period 1. Then X is a self- replicating multi-tiling 
family. 



Proof. Since the fundamental cycle oi X = (Xi,...,X„) has period 1, it contains a single 
element {xq}, and since it is irreducible, by Theorem |2.2| it has multiplicity one. Without 
loss of generality we assume that xq G Xi, so xq = Axq + d for some d G Dii. Let = 
[sf\...,S^n^) such that sf^ = {xq} and ah other sf^ = 0. Define 5^™) := ^(5^°)) = 
{S^^\. . . , 51""^). It follows from the expression for given in ( |6.1[ ) that 

n 

= \l (A'^'isf) + V^) = A™xo + V^. (7.1) 
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Suppose that (Ti,...,T„) is the set of self-affine multi-tiles corresponding to {A,Vij). Now 
by hypothesis Theorem |2.4| (i) holds, so by Theorem |2.4| (iv) each Tj has positive Lebesgue 
measure, and each Tj satisfies T° = Xi. We have 

n n 

(J (T,- + 5f = U (T, + A"^xo + Vf^). (7.2) 
It follows from K^{Ti) = U"=i(rj + V]{) that 

n 

A™(Ti+xo) = |J(r. + 5'f^ (7.3) 

The unions on the right side of (|7.3D are all measure- wise disjoint. Taking the limit m^oo we 
see that O := Uj=i(^j + ^j) is a packing of W^. 

It remains to prove that O is a tiling. The set Q. = IJj=i(^i + "'^i) is closed and satisfies 

n 

A{n) = y(A(T,)+A(X,)) 

n n 

j=l i=l 

n n 

i=l 3=1 

n 

= \J{T, + Xi) = n. 

i=l 

We now argue by contradiction, and suppose Q M'^. Since O is closed, there exists a ball B 
of radius e>0 such that Bnn = $, which yields (B) n A^ (Q) = A^{B)nn = 0. But A is 
expanding. So by taking N sufficiently large A^(i?) contains a ball of arbitrarily large radius. 
This ball is disjoint from $7, so it is not filled by any translate Tj +Xj, Xj € Xj and 1 < j < n. 
Therefore Xj cannot be a weak Delone set, a contradiction. Thus we have a tiling. I 



Remarks. (1) There exist irreducible self-replicating multi-tiling families whose fundamental 



cycles have period exceeeding 1, see Example \8.% Thus the condition of Theorem [7J| is only 
a sufficient condition (not a necessary condition) to give a self-replicating multi-tiling family. 

(2) Let X be an irreducible Delone set family which satisfies an infiation functional equation 
'(/'(A') = X with primitive subdivision matrix S has a fundamental cycle of period p. By 
Theorem 2^ and Theorem 3T the unique solution T := (Ti, ...,r„) consisting of compact sets 
of the associated multi-tile equation has Tj of positive measure, which are the closure of their 
interiors. The set T + X := ljr=i(^« + -^i) gives a partial g-packing of for some q < p, 
using the tiles Tj. That is, each point of M*^ is covered with multiplicity at most q off a set of 
measure zero, a set of infinite measure has multiplicity exactly q, and possibly another set of 
infinite measure has strictly smaller multiplicity. 

One expects that in many cases this construction produces ap-thick multi-tiling of R'^. We 
formulate the following problem: 
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Problem. Let X be an irreducible Delone set family which satisfies an inflation functional 
equation = X with primitive subdivision matrix S and is uniformly discrete and has a 

fundamental cycle of period p. Is it true that T + X := ljr=i(^« + -^i) always a g-thick 
multi-tiling for some q < pi If not, what are extra conditions needed to ensure it? 



8. Examples 

Example 8.1. (Substitution multiset with unbounded multiplicity function) Let A = [3] and 
= {0,1,2,3}, so that |Di,i| = 4 > 3 = |det(A)|. Take the seed 5^ = = {0}. Then 
the inclusion property S^^^ C S^^^ holds, hence 

X = Xi= hm 5f ^ 

defines a multiset X. The multiset X C Z>o and each point / € Z"*" occurs with finite 
multiplicity mx{l)- One can verify that mx(3'^) = k + 1, which shows that m{l) is unbounded 



as / oo. This example corresponds to case (ii) in Theorem 3.5 



Example 8.2. (Discrete substitution set family that is not uniformly discrete) Let A = [3] 
and Vi^i = {0,1, vr}, with vr = 3.14159.... Take the seed 5^°) = 5^^ = {0}. Then C ^^^^ 
and the limit 

Xi = lim 5f ^ 

exists. In this case the multiset Xi C M>o is discrete, and its elements all have multiplicity 
one. It is f 
all satisfy 



one. It is easy to show that it has linear growth. Indeed the 2-3" ^ elements in 5'|"''\S'|'^ 



3" < X < ^(3" + 3"-^ + • • • + 3 + 1) < 2^ • 3^^ . 
The associated multifile functional equation is 

A(ri) = Tiu(ri + i)u(ri + 7r) . 



The compact solution Ti to this equation has Lebesgue measure zero, see Kenyon [12| or 



Lagarias and Wang |17]. It follows from Theorem |6.1| that Xi cannot be uniformly discrete 



Example 8.3. (Inflation functional equation having infinitely many discrete solutions) Con- 
sider the inflation functional equation with A = [2] on M with n = 1 and Pi.i = {0, 1}, This 
satisfies the hypotheses of Lemma |3.7[ The associated tile is Ti = [0, 1]. The allowed starting 
points for a cycle Y of period p (an arbitrary positive integer) are given by ( ^.16 ), which gives 



m 

xo = for < m < 2P - 1. 

Each such xq generates an irreducible multiset X^^p consisting of a single multiset Xi = 
Xi{m,p) which has all multiplicites one. Now suppose that p > 2 and restrict to those 
Xo = — m/(2^ — 1) is in a cycle of minimal period p, with p > 2. This always happens when 
g.c.d.{m,2P — 1) = 1, and possibly in other cases as well, but rules out the cases m = 
and m = 2*^ — 1, which generate minimal cycles of period 1, having the symbolic dynamics 
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did2 ■ ■ ■ dp = 0^ and 1^, respectively. Thus — 1 < xq < 0, so that the tile Ti + xq includes in 
its interior, and the same holds for the other p—1 tiles in the periodic cycle. We conclude that 
the cycles of minimal period p generate a multiple tiling of M with thickness p, using copies of 
the tile Ti. 



Example 8.4. (Inflation functional equation having no nonempty discrete multiset solution) 
Consider the inflation functional equation with A = [2] in M, and n = 2 with digit sets 
Vii = Vi2 = 'D21 = V22 = {0, 1}. This data satisfies the hypotheses of Lemma 3^, but the 
inflation functional equation "il^{X) = X has no discrete multiset solution. If fact, all elements 
in X must have infinite multiplicity. To see this we iterate the inflation functional equation to 
obtain obtain 

X, = (Xi+P™)V(X2 + P-), i = l,2. 

Now each T>^ has an underlying set = {0, 1, . . . , 2™ — 1} with each element having multi- 
plicity 2™. Therefore the multiplicity of each element in Xj is at least 2"*. Hence no element 
in any Xi can have a finite multiplicity. 



Example 8.5. (Irreducible weak substitution Delone multiset family with bounded multiplici- 
ties, some multiplicities exceeding one) Consider the inflation functional equation A = [3] on 
M, with n = 2 and digit sets = {tt + 3}, X'1,2 = {!}, ^'2,1 = {-3,7r}, 1)2,2 = 0- The 

subdivision matrix S = \ \ is primitive, and its Perron eigenvalue A(S) = 2. We claim 

that the cycle y = {0 € Xi, 1 G X2} of period p = 2, both points of multiplicity 1, generates 
a discrete irreducible multiset family X = {Xi,X2). The only point reached in one step by 
exiting from the cycle is y = vr + 3 G Xi, which is reached in two different ways, one from 
G Xi and one from 1 G X2, so has multiplicity 2. All other points are descendents of y, and 
necessarily fall in the interval [tt + 3, 00). All the maps x — >■ 3x + d for d G "Dij are expanding 
outside the interval [—4,4] with expansion factor at least 1.4, so it follows that X = {Xi,X2) 
exists, and is irreducible and discrete. We claim that the descendents of y are distinct and 
consequently all have multiplicity two. To show they are distinct, note that all descendants 
of y have the form x = (3'^ao + 3^~^ai + ... + Sa^-i + afc)7r + r, in which ao = 1 and each 
subsequent = or 1, and r is an integer. The sequence (oq, ai..., afc) completely specifies 
the digit sequence leading to x, which determines which of Xi and X2 the point x belongs to; 
furthermore every such digit sequence is legal. Since all integers 3'^ao + 3^~^ai + ... + 3afc_i + afc 
are distinct, and vr is irrational, we conclude that all x are distinct. Thus all points in Xi and 
X2 have multiplicity two, except the generating cycle Y, whose two points have multiplicity 
one. 



Example 8.6. (Substitution Delone set that is not self-replicating) Let A = [3] and T>i^i = 
{ — 1,0, 1}. The associated tile Ti = [— ^, which consists of all balanced ternary expansions 

, d,- G {-1,0,1} . 



x = Y^ djS 

i=i 

The set 5^°) = {xq = -1/8, xi = -3/8} has 

5^°) C ^(5^°)) = {-17/8, -11/8, -9/8, -3/8, -1/8, 5/8}, 
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hence S^^^ generates an irreducible discrete multiset family X = {Xi} which consists of the 
single set Xi given by 

x,= Y^ s^"^ ■ 

k^oo 

The set Xi is irreducible and Sq is its generating cycle of period 2. A calculation gives 

x, = (J + z)u(-iuz) . 

It is a Delone set, and Xi + T is a multiple tiling of M of thickness 2. The thickness equals 
the period of the generating cycle, since both elements of S^^^ lie in the interior of the tile Ti. 
Thus Xi is a substitution Delone set but not a self-replicating Delone set. 



Example 8.7. (Self-replicating Delone set family having a primitive cycle of order larger 
than one) Let A = [—2] and Vn = {—2,-1}. Then X = (Xi) with Xi = Z is an irreducible 
substitution Delone set family whose fundamental cycle is {0,-1} and has period 2. The 
corresponding self-affine tile is Ti = [0, 1]. So T + X = Ti + Xi tiles M, hence is a self- 
replicating Delone set. 



Example 8.8. (Non-periodic and aperiodic self-replicating Delone set families) A Delone set 
X is (fully) periodic if it has a full rank lattice of periods Ax = {t : X = X -\- t}. It is 
non-periodic if it is not fully periodic and it is aperiodic if it has no periods, i.e. Ax = {0}. 
An example of a two-dimensional self-replicating Delone set, which not fully periodic, but 
has a one-dimensional lattice of periods, was given in Lagarias and Wang [18, Example 2.3]. 
Recently Lee and Moody |21| constructed many self-replicating Delone sets which are non- 
periodic, including aperiodic examples, whose points are contained in a lattice in W^. They 
give such examples associated to non-periodic tilings including the sphinx tiling of Godreche ||3| 
and the chair tiling. 
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